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^ ■ Abstract 

K^ ' We extend a recent billiard model of the nuclear A/^-body Hamiltonian to con- 

^— V I sider a finite two-body interaction. This permits a treatment of the Hamilto- 

f^ nian by a mean field theory, and also allows the possibility to model reactions 

a^ 

^^ I between nuclei. The density and the mean field potential can be accurately 

described by a scaling function which shows the qualitative features of the 
^ ' liquid drop picture of the nucleus. 

> 
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The nuclear many-body problem is often approximated by plienomenological mean field 
or liquid drop models which serve as good starting points to study the shape and the spec- 
trum of low lying single particle and collective excitations ||l]-§[ • This well-known approach 
is a consequence of the inadequate knowledge of the microscopic nuclear Hamiltonian, and 
the difficulty in dealing with the many-body problem. The mean field approximation is 
commonly employed to make the many-body problem tractable. Therefore it is of interest 
to question of how well mean-field theory works for a given microscopic Hamiltonian that 
can be solved exactly. 

A promising model for such a program may readily be obtained from a recent billiard 
model of the nuclear A^-body Hamiltonian. This is an interacting A^-body system with 
two-body interactions that is rich enough to show various features of self-bound many-body 
systems yet simply enough to allow for practical calculations and an understanding 0]. 

Generalizing the model of ref. |^, we consider the Hamiltonian 

N 2 

^ = E|r + En|r-;-r,|), (1) 

where r^ is a three-dimensional position vector of the i-th nucleon and Pi is its conjugate 
momentum and the interaction given by 



V{r) = < 



— Vn for r < a, 

(2) 
for r > a, 



where Vq is a positive constant. Total energy E, momentum P and angular momentum 
L are conserved quantities. Interactions alter the momenta of particles whenever they are 
distance a apart, and otherwise they move freely. 

For Vo = oo one obtains a billiard model, which was recently studied in ref. [^0. Due 
to the simple form of the two-body interaction, the numerical integration of the equations 
of motion is quite easy, scaling with particle number as A^logA^ rather than the usual 
A^^ typical of the classical A^-body problem. It was found that the dynamics is dominantly 
chaotic and ergodic. In the center of mass system one finds a constant single-particle density 



inside a circle of diameter a that quickly drops to zero within a thin surface region. The 
phase space structure does not depend on energy in such billiards. 

Finite values of Vq introduce an energy scale and a finite binding energy of the system. 
The Hamiltonian then permits reactions like particle capture or emission, giving it important 
features for a model of the nucleus and its dynamics. 

In this note we will examine the mean field theory of of the A^-body Hamiltonian (|l|,@). To 
this purpose we introduce the single-particle phase space density /(r, p). Let p^ /2'm + W{r) 
be the mean field Hamiltonian. As is well known, any density which is a functional of 
the mean field Hamiltonian solves the corresponding Vlasov equation 0. However, we are 
interested in the situation corresponding to thermal equilibrium and choose the canonical 
distribution function 

/(r.p1ocexp(-H!Z2!!lll!^), (3) 

where T denotes the temperature. Obviously, the phase space distribution function is a 
product of the normalized momentum space density 4(27r?TiA;T)^/^ exp (—p'^/2mkT) and 
the number density 

/ W(r)\ 
n(f) = cexpl ^j. (4) 

The normalization c is given by 

r r -1-1 

d^rn{r) 
The mean field potential ly (f^ depends on the number density and the two-body interaction 

W{r)= [ d^xV{\f-x\)n{x). (5) 



Eq.(H) together with eq.(|^) are the mean field equations to be solved self- consistently. 

We next observe that the solution of the mean field equations (§-^) will have an inter- 
esting scaling property. Let W{f) and nlf) denote the mean field potential and density for 
parameters N, Vq and kT. Then aW{r) is the mean field potential for parameters aNVo 



and akT, while an{r) is the density for parameters aN and Vo/akT. Thus, the normahzed 
quantities V{r)/NVo and n{r)/N depend only on the ratio 



7 = NVo/kT 



(6) 



of parameters. This scaling behavior comes directly from eqs. (^,^) and does not depend on 
the specific form of the two-body interaction (^. 

It is well known that the density of an open system can be represented as a thermal 
distribution function only for 7^1. In this case the density is exponentially small outside 
the nucleus and can be omitted. In what follows we restrict ourselves to this low temperature 
regime. 

We are interested in spherical symmetric solutions of the mean field equations. In this 
case eq.(BD may be written as 



W{r) = / dxn(x)v{r,x), 




(7) 



where 



v{r,x) = —27iVo < 





2x2 



for a < \r — x\ 
ioT a > r + X, 



(8) 



^ [a^ — (r — x)'^] for |r — a;| < a < r + x. 



2r 



For computational purposes it is useful to rewrite the mean field potential as 

W{r) = -27rVo 

+ e{a-r) [ dx^{ 
J 2r 

max(a— r,r) 

r -1 

+ 0(r-a/2) / dx—(a'^-(x-rY)n(x) 
J 2r 



a—r 

9(a — r) / dx2x'^n{x) 


a+r 



[a^ — (x — rV) n(x) 



\r-a\ 



(9) 



where 9(x) is the unit step function. 

We solve the mean field equations numerically by iteration. We take an initial density 
n{r) that is constant inside a sphere of diameter a and zero outside. Upon several iterations 



one obtains a converged solution. Fig. |l] shows the resulting number density n{r)/N and 
mean field potential W{r)/NVQ for different values of the parameter 7. For large 7 the mean 
field potential is fiat for r < a/2 and increases until r ^ 3a/2 while the density is constant 
for r < a and quickly drops to zero around r ^ a/2. Decreasing values of 7 (i.e. decreasing 
values of Vq and/or increasing temperatures T while the number of particles N is kept fixed) 
lead to a thicker surface region and to a more extended mean field potential. However, for 
r < a the potential is still a flat bottom potential and the density is roughly constant. Thus, 
the system exhibits qualitative features of the liquid drop model. However, our simple model 
does not yield a saturation of the density or the binding energy with increasing A^. This is 
due to the absence of a repulsive potential core. 

In the limit 7 -^ 00 the density approaches that of a sharp-edged liquid drop, constant 
inside a sphere of radius a/2 and zero outside. For finite values of 7, the surface thickness 
can be estimated assuming a trapezoidal density distribution. Taking the width parameter 
as e and demanding self-consistency in the neighborhood of r = a/2, one finds the relation 

e/a = (37)"^/^ (10) 

It is also interesting to compare the mean field results with the interacting system (0,0). 
Agreement may only be expected for sufficiently large numbers of particles Nq. In what 
follows we set A'"o = 300. The initial conditions of the A^o-body system are chosen such 
that the positions are drawn at random inside a sphere of diameter a while the momenta 
are drawn from a Maxwell-Boltzmann distribution with temperature T. We follow the time 
evolution of the many-body system for roughly 10^ interactions to allow for equilibration. 
At the end of the evolution we take the positions in the center of mass system and compute 
the integrated density 

r 

7V(r)=4./dxx^„(x) (11) 



by counting the number of particles inside a sphere of diameter r. Fig. |^ compares the 
average of ten runs with the mean field result for various values of the temperature. The 
agreement is rather good confirming the validity of the mean field approximation. 



Finally, we ask what values the parameters should have to correspond to physical prop- 
erties of nuclei. Let us consider a nucleus of mass A ^ 75. Then the radius is given by 
a/2 ^ 1.2 A^^^ Ri 5 fm. The liquid-drop surface thickness is roughly e ^ 1 fm, yielding 
from eq. ( [T0| ) 7 ~ 30. We would like to set the energy scale by the depth of the mean 
field potential, which should be A^Vq ~ 50 MeV to correspond to a typical Woods-Saxon 
potential of nuclear physics. Turning to eq. (7), we see that the two parameter values can be 
obtained taking the temperature as kT ^ 2 MeV. Equilibrated nuclei at such a temperature 
can be easily produced in heavy ion reactions, so the model might have some applicability 
to nuclear reactions. However, in detail the Fermionic nature of the nuclear many-body 
problem will make considerable differences from the present classical model. In particular, 
the chemical potentials at a given temperature are very different, which would be important 
for the nucleon evaporation rates. 
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FIGURES 




FIG. 1. Mean field potentials (lower panels) and densities (upper panels) for various values of 
the parameter 7 = NVo/kT. 
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FIG. 2. Integrated number density from microscopic calculations (dashed lines) compared to 
mean field results (solid lines) for different values of the parameter 7 = NVo/kT 



